Non-perturbative predictions for cold atom Bose gases with tunable interactions 
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We derive a theoretical description for dilute Bose gases as a loop expansion in terms of composite- 
held propagators by rewriting the Lagrangian in terms of auxiliary helds related to the normal and 
anomalous densities. We demonstrate that already in leading order this non-perturbative approach 
describes a large interval of coupling-constant values, satisfies Goldstone's theorem, yields a Bose- 
Einstein transition that is second-order, and is consistent with the critical temperature predicted in 
the weak-coupling limit by the next-to-leading order large-N expansion. 
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Nearly a century after the first observation of the 
lambda transition in liquid hclium[l], a quantitative, 
first-principles description of strongly-correlated bosons 
remains a challenge. After the transition was recog- 
nized as the onset of superfluidity [2], the connection with 
Bose- Einstein condensation (BEC) was proposed [3], but 
it was Bogoliubov's work[4] pointing out that the disper- 
sion of the elementary BEC excitations satisfy the Lan- 
dau criterion for superfluidity[5] that motivated weakly- 
interacting BEC studies to investigate superfluid proper- 
ties. In weakly-interacting systems, the many-body prop- 
erties do not depend on the shape of the interaction po- 
tential, but only on the s-wave scattering length, do, and 
the boson fluid acts as point-like interacting particles [6]. 

Unlike liquid helium, cold atoms remain point-like even 
when the scattering length is tuned near a Feshbach reso- 
nance. Then, strongly-correlated cold atom bosons offer 
the exciting prospect of studying point-like strongly in- 
teracting bosons, possibly in the universal regime where 
the scattering length greatly exceeds the inter-particle 
distance and the latter becomes the only relevant length 
scale [7] . This hope appeared thwarted when it was shown 
that the three-body loss rate in cold atom traps scales as 
o,q near a Feshbach resonance[8]. In accordance, the uni- 
versal regime was reached only in ultra-cold fermionic 
gases [9], where the three-body loss is reduced by virtue 
of the Pauli exclusion principle. However, the recent ob- 
servation that three-body losses are strongly suppressed 
in optical lattices when the average number of bosons per 
site is two or less [10], rekindles the prospect of studying 
medium and strongly-correlated cold atom bosons. Novel 
cold-atom trap technologies that produce stable, flat po- 
tentials bound by a sharp edge[ll], suggest the study of 
finite-temperature properties such as the BEC transition 
temperature T c and the superfluid to normal fluid ratio 
and depletion, at fixed density, p. 

At finite temperature, the description of BEC's re- 
mains a challenge even in the weakly-interacting regime. 
Standard approximations such as the Hartree-Fock- 
Bogoliubov and the Popov schemes, generally fall within 
the Hohenberg and Martin classification[12] of conserv- 



ing and gapless approximations, which implies that they 
either violate Goldstone's theorem or general conserva- 
tion laws[13]. These approximations generally predict 
the BEC transition to be a first-order transition, whereas 
we expect the transition to be second order [14]. 

In this paper, we present a new theoretical framework 
that describes a large interval of p 1 / 3 ao-values, satisfies 
Goldstone's theorem and yields a Bose-Einstein transi- 
tion that is second-order, while also predicting reasonable 
values for the depletion. Furthermore, this framework 
can predict all experimentally relevant quantities within 
the same calculation, determining fully consistently 
quantities such as T c , the collective mode frequencies [15] 
and the compressibility (which characterizes the density 
profile in a shallow trap[16]). In contrast with other re- 
summation schemes, such as the large-iV expansion[17] 
or functional renormalization techiques[18], here we treat 
the normal and anomalous densities on equal footing. 

In our approach, we generate a one-parameter family 
of equivalent Lagrangians. We choose this parameter to 
reproduce the one-loop result at mean-field level in the 
weakly-interacting limit. Thus, we identify the optimal 
auxiliary-field Lagrangian for the purpose of a systematic 
non-perturbative expansion. Then, the critical temper- 
ature variation in leading order is the same as the one 
found in the next-to-leading order large- N expansion. 

In dilute bosonic gas systems, the classical action is 
given by S[<f>, (fr* ] = Jdx C[ 4>, 4>* ] , with dx = dt d 3 x and 
the Lagrangian density 

c[<f>,f] = * [<W & 000) - (d t <p* (*)) 4>(x) } 

- <t>*(x) { - m} <t>(x) - \ I 0W| 4 . (1) 

Here, /i is the chemical potential and the coupling is 
A = Anh 2 ao/m. To account for the contributions of the 
normal and anomalous densities, we use the Hubbard- 
Stratonovitch transformation[19] to introduce the real 
and complex auxiliary fields (AF), x( x ) an( i A{x). We 
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add to Eq. (1) the AF Lagrangian density[20, 21] 



£ aux [<M*,X,^*] = ^ [x(x)~\cosh6\<f>(x)f 



2A 



A(x) - A sinh6^ 2 (a;) 



(2) 



where 9 is the mixing parameter between the normal and 
anomalous densities, x( x ) ancl A(x). The usual large-N 
approximation [21] is obtained when 9 = 0. Then, the 
action becomes 



S[$,J} = S[4> a ,x,A,A*,j a ,s,S} 
1 

~ ~2 



(3) 



dx dx' </> a (x) G- la b [ X , A](x, x 1 ) 4>\x') 

+ J dx{[ X 2 (x)-\A(x)\ 2 ]/(2X)- S (x) X (x) 
+ S*(x)A(x) + S(x)A*(x) + j*(x)</>{x) + j(x)<f>*(x) } , 
with 

G- la b [ x ,A] = {Gu la b + V a b [x,A}(x)}6(x,x') , 
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V\[ X ,A](x) = 



X(x)cosh8 ~A(a;)sinh6 
-A*{x) sinh 6* X { x )cosh.9 



Here, we introduced a two-component notation with 
c/) a (x) = {(/>(x),4>*(x)} for a = 1,2. $(x) and J{x) sig- 
nify the five- component fields and currents. The gener- 
ating functional for connected graphs is 



Z[J] = e lW W h =AT J D$e 



iS[<S>;J]/h 



with S[<&; J] given by Eq. (3). Performing the path inte- 
gral over the fields (f> a , we obtain the effective action 

eS cS [ X ; J,e] = \U dxdx'j a (x)G[x] a b (x,x')j a (x) 



dx ■ 



Xi(x) x l (x) 
2A 



Si{x) X l (x) 



2i 



-Tr 



HG- 1 }}, 



where X l {x) = {x(x),A{x)/V2,A*(x)/V2}, S l (x) = 
{s{x),S{x)/V2,S*(x)/V2}. The small parameter e al- 
lows us to perform the remaining path integral over x % 
using the stationary-phase approximation. As shown in 
Ref.20, e counts loops in the AF propagator in analogy 
with h, and provides the loop expansion of the effective 
action in terms of x propagators. Next, we expand the ef- 
fective action about the stationary points, Xo( x ), defined 
by 5S c e[x\ j]/ 5xi{x) = 0. Hence, we obtain 



Xo{x) 
A 

A (x) 



= {\M X )\ 2 + ^.^[G(x,x)}} cosh0 + s(x) 



x {<j) 2 (x) + -.G 2 !(x,x)} smhO + S(x), 



where we introduced the notations 

r [Xo](x) = J dx'G[ X o] a 6 (^')jV)- 

We emphasize that both xo and Ao include self-consistent 
fluctuations. Expanding the effective action about the 
stationary point, we write 

s cS [x; J] = SMxo;J} + l JJ dWa/i^xo] 

x[ X i (x)-x (x)]{x j (x , )-xi(x , )]+--- , (5) 
where D~^{ X ^ x 1 ) is given by the second-order derivatives, 

s 2 s cS [x a ] 



D^[xo](x,x') 



5x i (x)Sx j (x / ) 



evaluated at the stationary points. By keeping the gaus- 
sian fluctuations and Legendre transforming, the one- 
particle irreducible (1-PI) graphs generating functional 
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r[$] = / dxj a {x)4> a {x)-w[J] 

' 1 f dxdx' <t> a (x)G- 1 \x] a b(x,x')<f> b (x') 



Xi{x)x l {x) fi T/| r „_ lf s , 
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is the negative of the classical action plus self-consistent 
one-loop corrections in the <j> a and Xi propagators. 

To leading order in the AF loop expansion (LOAF), 
one sets e = in the right-hand-side of (6). The static 
part of the effective action per unit volume is 



V cS m = (xcoshfl _ ) \cf)\ 2 - 1 (A* <f> 2 + A0* 2 ) sinht 



x 2 -W 



■^{HG-'ix]}} 



(7) 
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Translating (7) to the imaginary time formalism, we find 

^H<r'Wl-/^{f + Jm-^]}, 

where u: 2 . = (e^ + x cosh6> — ^) 2 — \ A\ 2 sinh 2 9 and tk = 
k 2 /{2m). At the minimum, we have 



6V eS [$] 



= (xcashd- (ji)<I> - A Biahefa =0. (8) 

Using the U(l) gauge symmetry, we choose <fto to be 
real. Then, A is real and the dispersion, ui 2 = 6k(ek + 
2A sinh 9), represents the Goldstone theorem. Next, we 
set sinh# = 1, such that LJf. reduces to the Bogoliubov 
dispersion, u>k = \J £fe(e/c + 2A 0q), in the limit of vanish- 
ing quantum fluctuations in the anomalous density. We 
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FIG. 1. (Color online) Normal density, x\ an< i anoma- 
lous density, A, from the LOAF and PA approximations, for 
p 1//3 ao = 1. T c and T* indicate vanishing condensate den- 
sity, po, and anomalous density, A, respectively. PA leads 
to a first-order phase transition, whereas LOAF predicts a 
second-order phase transition. We have that T c — T* in the 
PA, but not in LOAF. In LOAF x' and A are equal until T c . 



note that the leading-order (LO) in the large-N expansion 
corresponds to 9 = 0. This leads to the noninteracting 
(NI) dispersion, ui k = e k , and we conclude that the large- 
N expansion is not a suitable starting point, because it 
is incompatible with the Bogoliubov spectrum. 

Using standard regularization techniques [22], the 
renormalized effective potential is written as 



d 3 k 
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w fe - e k 



(A*4> 2 + A<f> 
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where x' — V%X ~ A 4 an d °~ 
The gap equations, obtained from 514ff[3>]/<5x 

d 3 fc r 1 + 2n(w fc ) 
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where n(tu k ) = [exp^^/Ac^T) — is the Bose- Einstein 
particle distribution. At the minimum of the effec- 
tive potential we have, (x'o ~ ^o) 0o = 0, see Eq. (8), 
and we replace p by the physical density using p = 
-dV e &l$o]/d(i = (xo + M)/( 2 A). The density is used to 
rescale Eqs. (9), and the ensuing phase diagram problem 
depends only on the dimensionless parameter, p x / 3 ao, 
and the coupling constant becomes A = 8n p 1 / 3 ^. In the 
broken symmetry phase, we have x'o = and the dis- 
persion relation, cj| = e k (e k + 2%q). The condensate 
density is denoted by po = <pQ. At weak coupling and 
T = 0, our results coincide with the Bogoliubov (one- 
loop) approximation [14], p — 8tt pa [l + (32/3) \/ pa^j-K ] . 
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FIG. 2. (Color online) Temperature dependence of the con- 
densate fractions from LOAF and PA, compared with the NI 
result, for p 1//3 ao = 0.1 and p 1 ^ i aa = 1. Because at T c the PA 
and NI dispersion relations are the same, PA does not change 
T c relative to the NI case. LOAF increases T c . 



0.40 



h? 0.30 



< 



0.15 




P" J a„ 



FIG. 3. (Color online) Relative change in T c with respect to 
NI, as predicted by LOAF as a function of p 1//3 ao. The inset 
shows that in the weak-coupling regime, LOAF produces the 
same slope as the next-to-leading order large-N expansion [17]. 



We compare the LOAF results with the predictions of 
the Popov bosonic approximation (PA) [23] . PA is gener- 
ally recognized as an accurate theoretical description of 
experimental data in weakly-coupled dilute trapped Bose 
gases [24], as long as the densities of the condensed and 
noncondensed atoms are comparable with each other. 
Unfortunately, PA produces an artificial first-order phase 
transition at T c . Formally, PA is obtained from Eq. (9) by 
setting Aq = Xo — ^/°o an( i neglecting the quantum fluc- 
tuations in the anomalous density. With this substitu- 
tion, the PA dispersion relation reads w 2 , = e k (e k +2Xpo). 

In Fig. 1 we depict the temperature dependence of the 
normal density x' i arL d anomalous density, A, at con- 
stant p 1 / 3 ao, as derived using the LOAF and PA approx- 
imations. For illustrative purposes, we set P 1/3 a = 1 
and the temperature is scaled by its NI critical value, 
T = (27rft 2 /m)[p/C(3/2)] 2 / 3 , where ((x) is the Riemann 
zeta function. We identify two special temperatures, at 
T c where the condensate density vanishes, and at T* 
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where the anomalous density, A, vanishes. These tem- 
peratures are the same in the PA formalism, but they 
are different in LOAF. The existence of a temperature 
range, T c < T < T* , for which the anomalous density, 
A, is nonzero despite a zero condensate fraction, <j>, is 
a fundamental prediction of LOAF. In this temperature 
range, the dispersion relation is expected to depart from 
the quadratic form predicted by the Popov approxima- 
tion for T > T c . Above T c the solution of the PA equa- 
tions becomes multivalued, indicating that the system 
undergoes a first-order phase transition at T c . In con- 
trast, LOAF predicts a second-order transition. 

The temperature dependence of the condensate frac- 
tion, po/p, is depicted in Fig. 2 for two constant val- 
ues of the dimensionless parameter p 1 / 3 ao, together with 
the NI result, p Q /p = 1 - (T/T ) 3 / 2 . Again, we observe 
that LOAF exhibits the correct second-order BEC phase 
transition behavior. Moreover, PA does not change T c 
relative to the NI case, because in the PA case we have 
T c = T* and the PA and NI dispersion relations are the 
same at T c . The LOAF approximation predicts an in- 
crease of T c compared with the NI case. 

As illustrated in Fig. 2, the LOAF and PA predictions 
may differ greatly even for temperatures, T <C T c . These 
differences are enhanced by a strengthening of the inter- 
action between particles in the Bose gas (a larger value 
ofp^ao indicates stronger coupling). The leading-order 
AF formalism produces a more realistic set of observ- 
ables away from the weak-coupling limit because of its 
non-perturbative character. In contrast, PA is appropri- 
ate only in the case of a weakly-interacting gas of bosons. 
The former is made explicit by studying the LOAF pre- 
diction for the relative change in T c with respect to To, 
as a function of p 1 / 3 ^. The inset in Fig. 3 demon- 
strates that in the weak-coupling regime, / o 1 / 3 ao <C 1, 
LOAF produces the same slope of the linear departure 
derived by Baym et al. [17] using the large-N expansion, 
but at next-to-leading order. The LOAF corrections to 
the critical temperature are due to the inclusion of self- 
consistent fluctuations effects in the mean-held x' an d A 
densities. A summary of AT c /Tq theoretical predictions 
is found in Ref.14. For p 1 ^ 3 a 3> 1, LOAF predicts that 
AT c /To — > 0.396 when the system approaches the unitar- 
ity limit. Despite that most current experiments probe 
only the p 1 / 3 ao <?C 1 regime, future experiments [11] may 
access the medium-to-strongly interacting regime, and 
verify this non-perturbative prediction. 

One can systematically improve upon the LOAF ap- 
proximation by calculating the 1-PI action order-by- 
order in e. The broken U(l) symmetry Ward identities 
guarantee Goldstone's theorem order by order in e [20]. 
For time-dependent problems, however, this expansion is 
secular[25], and a further resummation is required. The 
latter is performed using the two-particle irreducible (2- 
PI) formalism[26]. A practical implementation of this ap- 
proach is the bare- vertex approximation (BVA) [27] . The 



BVA is an energy-momentum and particle-number con- 
serving truncation of the Schwinger-Dyson infinite hier- 
archy of equations obtained by ignoring the derivatives 
of the self-energy, similarly to the Migdal's theorem[28] 
approach in condensed matter physics. The BVA proved 
effective in the case of classical and quantum A</> 4 field 
theory problems[29] and can be applied to the BEC case. 

To summarize, in this paper we introduce a new non- 
perturbative resummation formulation for the BEC prob- 
lem. At mean-field level, this approach meets three im- 
portant criteria for a satisfactory mean-field theory for 
weakly-interacting bosons[14]: i) the excitation spectrum 
is gapless (to preserve Goldstone's theorem), ii) LOAF 
reduces to the known results from Bogoliubov theory at 
T = and weak coupling, and iii) predicts a second-order 
BEC phase transition. The latter suggests that a AF for- 
mulation of the Lagrangian for systems of cold fermionic 
atoms may also impact the study of the BEC to BCS 
crossover in dilute fermionic atom systems [30]. 

Work performed in part under the auspices of the U.S. 
Department of Energy. The authors would like to thank 
E. Mottola and P.B. Littlewood for useful discussions. 
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